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Abstract 
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by Buscher's T-duality transformation. The computation of the two loop 
perturbative free energy density is performed in the case of a certain defor- 
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sional regularization and the geometric sigma model perturbation theory. 
We obtain agreement between the free energy density expressions of the 
two models. 
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1 Introduction 



Among the wealth of different dualities relating the perturbative string theories 
and M-theory, the most 'ancient' one is T-duality. At the same time T-duality 
was the starting point in the discovery of D-branes [Q]. In string theory T- 
duality can be proven in arbitrary order of the string perturbation theory |3|, 
as long as the vacuum preserves conformal invariance @. Generalizing the 
R — > symmetry of the toroidally compactified strings, Buscher |p gave a set 
of equations that describe the transformation of the Neveu-Schwarz background 
fields of the sigma model action, and works whenever there is an isometry. 
Denoting by gfj_p and bfj_i, the background metric and the antisymmetric tensor 
field, the explicit form of the Buscher transformation is (since we work on a flat 
world sheet, we neglect the dilaton): 

900 — ) 90a — , Ooa — , 

500 900 900 

90a90/3 — boabof3 7 , 50a^0/3 " 50/3^00 
9al3 — 9al3 , 0^(3 — Oa/S , 

900 500 

where g^u and b^u denote the new background fields. 

To derive this transformation Buscher used functional integral arguments, 
that in the meantime have become widely known, and applied in many context 
(see e.g. [^). The idea is to gauge the aforementioned isometry of the sigma 
model action and impose a constraint using a Lagrange multiplier. Integrating 
the multiplier one recovers the original theory to start with, while after a gauge 
fixing and integrating over one of the original fields gives the dual theory. 

It was also shown [^, that, at the classical level, the duality transformation 
rule can be recovered in an elegant way by performing a canonical transforma- 
tion ||8|. This clearly shows that the models connected by the Buscher trans- 
formations are equivalent classically. At the quantum level the only loophole 
of the path integral argument mentioned above is that it neglects the effect of 
the regularization and renormalization. Though it was shown in |Q] that for 
conformal invariant models one has full quantum equivalence, in other words 
the path integral argument holds, this is not the case in the non-conformal 
setting. It was shown in |^] for a deformed SU (2) principal sigma model, and 
further clarified in |jl^ for several different deformed SU{3) principal sigma 
models, that Buscher's formula - as applied to renormalized quantities - has to 
be modified to give two loop quantum equivalence. 

Later it was shown that the one loop beta functions of the original and dual 
models always agree [11 1. Work has been done to establish the corrections to 



Buscher's formulae |12, 13, 14 1. Advances were made from a different point 
of view. It is widely known that the low energy degrees of freedom of the 
sigma model can be described using an effective action, that contains gravity 
in target space. This fact constrains the possible leading terms of the low 



energy effective action to a computable form, that is known. It was shown [15| 
that the low energy effective action consistent with the two loop sigma model 
beta function equations is not invariant under Buscher's transformation. The 
leading part of the above action, the one loop part, nevertheless is invariant, 
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in accordance with the one loop findings of [|lT[. Using the non-invariance of 
the two loop action it was found how to modify Buscher's transformation, with 
order a' terms, such that the two loop action would remain invariant under the 
modified transformations. What is not clear after all is how to pull back the 
modification of the Buscher's formulae found in the low energy action to the 
sigma model. Related work in the supersymmetric case, concentrating on the 
absence of the above mentioned corrections, was done in |T^. 

Any satisfactory criterion of the quantum equivalence among dually related 
sigma models should be based on the comparison of physical quantities as op- 
posed to just considering beta functions. If there are global symmetries in the 
model then their associated conserved quantities (Noether currents) may be 
considered physical. The definition of physical quantities, however, is not very 
clear in diffeomorphic invariant sigma models without a sufficient number of 
isometrics. To circumvent this problem the study of Weyl anomaly coefficients 
was suggested in [^]. In the present paper - as an alternative - we study a 
thermodynamic quantity namely the free energy density in the presence of a 
chemical potential in the dually related sigma models. This quantity surely 
qualifies as physical, thus its equality in the two models gives a non trivial 
check on their quantum equivalence. Furthermore the free energy density can 
be computed perturbatively - at least in asymptotically free models, thus one 
can compare the two free energy densities using the first few orders of pertur- 
bation theory. The aim of this paper is to carry out this comparison in the two 
loop order, where the first really 'quantum' effects appear, thus improving the 



almost 'classical' one loop case studied earlier in [|16|], [|17]] . 

The paper is organized as follows: in Section 2 we give a brief review of the 
pertinent facts that we need from the renormalization of the SU{2) principal 
a model and its dual and develop a Lagrangian with more parameters, that 
can accommodate both Lagrangians as special cases. In Section 3 we investi- 
gate the conserved charges that can be coupled to the models, calculate the 
Hamiltonian, find the ground states in the presence of the external field, and 
make the Lagrangians suitable for a perturbative computation. Section 4 deals 
with the definition of the perturbative free energy density, and its computation, 
using dimensional regularization. First we define the perturbative free energy 
density, and set up a scheme to compute it systematically to any order. The 
rest of Section 4 deals with the actual computation of the bare free energy up 
to two loops. In Section 5 we deal with the issue of renormalization and obtain 
the one and two loop renormalized free energies of the original and dual models, 
compare them, and improve them using the renormalization group. We close 
Section 5 with the analysis of the composite operator renormalization of the 
relevant operators. We make our conclusions in Section 6. 

2 Lagrangians and T-duality 

2.1 The deformed SU{2) principal a model 

This section has a twofold role. Primarily it is intended to give an overview of 
the results that we need in the rest of the paper and at the same time fix the 
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notations. Secondarily it extends some of the earlier results in a way suitable 
for our applications. 

In the following one parameter deformation of the SU{2) principal a- 
model was considered : 
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^o = -^{T.Jp''''+94j'''), (2) 



a=l 



where = G~^d^G = J^t"', and t'^ = a°-/2 with a"" being the standard 
Pauli matrices. Thus G is an element of SU{2) and g is the parameter of the 
deformation. From the Lagrangian @) it is clear that the global SU{2)l x 
SU{2)ji symmetry of the undeformed principal cj-model is broken to SU{2)l x 
U{1)r by the Jf^J^^ term. Setting g = Q corresponds to the principal cr-model, 
while for g = —1 the 0(3) cj-model is obtained as it can be seen from eq. (Q) 
below. 

The authors of 1^] investigated the renormalization of A and g in the two- 
loop order of perturbation theory, treating A as the coupling constant and g as 
a parameter. Using the Euler angles {(f),9,tl') to parameterize the elements of 
SU{2), G is written as 

Using this converts the Lagrangian of the deformed a model, which for the time 
being we shall call 'the original model', into the following form: 

+ (1 + g){d^ijf + 2(1 + g)d^m cos 6} . 

Using the Killing vectors of the SU{2)]^ x C/(1)r symmetry and exploiting 
the manifest target space covariance of the background field method it was 
proved in that the model is renormalizable in the ordinary sense: there is 
no wave function renormalization for 9, (j) and -0, while the coupling constant 
and the parameter got renormalized according to: 

Ao = /u~'Za(A,5)A, go = Zg{X,g)g. (5) 

Both in the one and in the two loop orders the residues of the single poles in 
ZxiX,g) = l-yx{X,g)/e+. . . and Zg{X,g) = l-yg{X, g)/e+. . . were determined: 

Note the sign difference between our formulas (^) and (^ , and the corresponding 
ones in . It is consequence of the fact that in our notation notation n = 2 + e 
rather than n = 2 — e as used in [Btl. 
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The standard definition of the (3 functions: /3q = A*^, l^-^j = /^s^> lead to 
the following two-loop [3 functions (eq. (20) in Ip): 

I3,=^a(i + g){i + ^{i-a))^ 

It is easy to see, that the 5 = resp. the g = —1 lines are fixed lines under 
the renormalization group flow, and (3\ reduces to the [3 function of the principal 
c7-model, resp. of the 0(3) cr-model on them. In the (A > 0, g < 0) quarter of 
the (A, (7) plane the renormalization group trajectories run into A = 0,5 = — 1; 
while for > they run to infinity. This implies that the g' = fixed line 
corresponding to the principal cr-model is 'unstable' under the deformation. 

The Lagrangian of the deformed cr-model, eq. (Q), exhibits two obvious 
Abelian isometries that can be used to construct two different (Abelian) duals: 
namely the translations in the (f) and il) fields. We call the models obtained this 
way the '(/> dual' respectively the 'V' dual' of the deformed a model (^. 

In Q it was found that for the '"0 dual' model, as summarized below, the 
renormalization of the coupling and the parameter are equivalent to that of 
the original model. Therefore, in the present context we deal with the original, 
deformed SU{2) principal cr-model and it's '^/j dual', which we shall simply call 
'the dual model'. 



2.2 The 'ijj dual' model 

For the Lagrangian of the -dual' model, using Buscher's formulae |5|, |^ , one 
has an expression analogous to eq. (^): 

Cd =i {{^^,9f + (d^^f sin2 6 + {d^x? + 2a cos Oe^^'d^xdu^) ■ (8) 

Here x denotes the (appropriately scaled) variable dual to ip, and (A,^) stands 
for the couplings of the dual model. One can show that Cd exhibits the expected 
SU{2) X U{1) symmetry for all values of the parameter a. 

The couplings of the original, (Q), and of the dual models, (^), are related 
(at the classical level) as a direct consequence of T-duality as 

A = A , a = x/l + g , g = g. (9) 

These relations can be maintained at the two-loop level if one performs the 
renormalization and, in addition, in both theories the couplings {g , g) are ex- 
pressed in terms of the corresponding renormalization group invariant quanti- 
ties, which at the end are set equal to each other In addition one also has 
to take into account the freedom in the choice of the renormalization group 
invariant, due to the scheme dependence of the two-loop /? function. 

The two-loop renormalization invariants, that characterize the fiows under 
the corresponding sets of /5 functions, for the original and dual models can be 
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easily computed: 



*'/f ;v (10) 

The next step is to express g (resp. a) in terms of the renormahzation group 
invariant. This can be done by inverting eq. ([To|) perturbatively. In both cases 
at leading order one has to solve quadratic equations and the next to leading 
order terms correct the results including the two- loop effects. Assuming that 
neither Mq nor Md vanishes, one obtains (the leading terms were obtained in 
Q, but we shall need the whole expression): 

The sign ambiguity can be removed by studying the renormahzation group 
flows, as it was briefly mentioned in connection with eq. (0). In the original 
model, it turns out that the interesting region is the vicinity of g = —1 with 
g > —1. In the dual model, since a = y/1 + g, choosing g > —1 one is uniquely 
led to the solution with the plus sign. Thus in both cases one has to consider 
the solution with the plus sign. Moreover if one sets Mq and M^) equal, Mq = 
Md = M, then the classical a = ^/l + g relation can be maintained. We note 
(though we didn't display it in eq. (pA])) that the two expressions differ already 
at the order of A^, as expected, since M is renormahzation group invariant only 
up to two loop order. 

Using g{\,Mo) {a'^{X, Mij)) in the two loop beta functions of the coupling 
constants of the original and dual models yields a universal expression |^]: 

As far as the coupling constant renormahzation is concerned, this universal beta 
function shows that the two models are equivalent, both are asymptotically free, 
and the actual value of M effects only the two loop coefficient. 

According to eq. ( pT|) one can express a in terms of the common renormal- 
ization group (RG) invariant. For the sake of simplicity it is useful to introduce 
the following notations: 

A = VX , ai = -r7= , ao = — 71= ■ (13) 

It will turn out that A = -v/A is the proper coupling of the two models. In terms 
of these 

a = aiA - a2A^ + o{A^) , g = -1 + ajA^ - ^afA^ + oiA^) . (14) 
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2.3 Unified description 



In order to make the computations more general it is useful to express both 
theories as particular cases of a generalized a model. To achieve this we rescale 
the field "0 — > ^i+g ^ ^^'^ introduce a = + g in place of g in the original 
model. After some obvious changes of symbols, Co and Cd can be described 
by the following unified a model: 

C =^{{df,ef + {df,cl))\l + r cos2 6) + (a^V)' + 2aa;^, cos 0d^cl)d''iP} . (15) 

The previous models can be recovered by the special choices: r = g and uif^u = 
r]^y for the original model, and r = — 1 and ujfj_y = — e^jy for the dual model. 
Expressing r in terms of the RG invariant parameter we have: 

r = -l + 2/?A2+ 7^^ + 0(^*5), (16) 

where (3 = a\/2 and 7 = — ^af in the original resp. /3 = 7 = in the dual 
model. Observe that both f3 and 7 are renormalization group (RG) invariant. 

This unified Lagrangian (^) is more than it might appear at first sight. 
As it was shown in Q and argued above, the deformed principal a model 
and its dual can be viewed as being quantum equivalent from the point of 
view of the two loop beta functions. The unifying Lagrangian (|l^) can be 
viewed as a genuine quantum generalization of the deformed principal a model 
and its dual. It reduces to the latter ones at special values of the parameters 
/3, 7 and ujfj_y, and at the same time gives the corresponding beta functions. 
Thus the renormalization properties of the two models are encompassed in this 
generalized Lagrangian. Based on this we will be able to renormalize the free 
energy in both theories at the same time, shortening the computations and 
obtaining better control on the different contributions.. Checking the quantum 
equivalence will amount to compare the renormalized free energies computed 
at the special values of the parameters (3, 7 and oj^u corresponding to the two 
models. 



3 The ground state and perturbative Lagrangian 

3.1 Outfine of the method - Noether currents 

So far testing the quantum equivalence of the dual models was mostly reduced 



to comparisons of the corresponding beta-functions 10, 11 1. Of course there 
are more to test before one can ascertain about an equivalence. In this paper to 
test the physical equivalence between the original and dual theories we couple 
both of them to some particular conserved charge Q. This is accomplished 
by modifying the respective Hamiltonians Hq [Hd) to Hq — hQ {Hq — hQ), 
where h is an external field (chemical potential type of parameter), having 
mass dimension one. The corresponding changes in the ground state energy 
densities (i.e. in the free energy densities, that we shall call for simplicity free 
energies) can be computed, at least in principle, to any order in perturbation 
theory. The comparison of the free energy densities, as functions of h, up to 
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a certain order of perturbation theory, in the original and dual models, then 
provides a useful check whether the two models do really, physically correspond 
to each other. For a comparison of this type to make sense both theories must 
be asymptotically free (to guarantee that perturbation theory applies), and of 
course we have to choose Q to be really the same. 

As the global symmetry group of both the original and the dual models is 
SU{2) X U{1) we may think naively that any linear combination of an SU{2) 
Noether charge and the U{1) Noether charge, x"Q° + yQu(i)-, can be used as 
the charge above Q to couple to the Hamiltonians. However for our test we 
need the same Q coupled to Hq and Hd, thus we can choose only such charges 
that are mapped to themselves under the canonical transformation connecting 
the original and dual models. This, of course, implies that the charge Q must 
stay local under the canonical transformation. 

We first point out the relation between the appropriate Noether currents 
and charges of the global symmetries of Cq and Cd- An exhaustive treatment 
will be given elsewhere [20|. Let's start remarking that under Abelian duality 
transformations the image of the U{1) current, that belongs to the distinguished 
isometry used in duality, is a topological current built from the dual field. Thus 
the image of the U{\)pi current of the deformed sigma model is the topological 
current of the x held, e^^d^X^ ^-iid as such, its charge should vanish on a topo- 
logically trivial 2d space-time. Therefore Qu{i) cannot be present in both Hq 
and H£). 

The next simplest possibility is to use the U{1) charges corresponding to 
the (/)-translations in both the original and dual theories: Nq and Nq. It can be 
shown (see [20| for details) that the canonical transformation which implements 
Abelian duality, effectively exchanges only and x' iPx ^'^^ "0') and leaves p^, 
pe unchanged. It is obvious that Nq is 'identical' under duality to Nq, but it 
is not entirely trivial (though it is true) that the space component of A'^^ really 
becomes the space component of iV^. 



3.2 From the Hamiltonian to the perturbative Lagrangian 

In conclusion, to test the physical equivalence of the original and dual mod- 
els, in both of them we introduce the coupled Hamiltonian densities 7^o,D = 
TCo,D — ^A^Q .[| Performing the inverse Legendre transformation on these quanti- 
ties one obtains the Co,d Lagrange densities of the coupled models. By explicit 
computation one can show that this procedure of obtaining Co,D is equivalent 
to the following formal gauging {d^C — > -D^C*) of the original Lagrangian Co,D' 

D^^p = d^ip, D^cp = d^^ + h6o^, D^e = d^e. (i7) 

This gauging, of course, can be done on the universal, common Lagrangian 
(P^); the outcome being denoted by C. As it can be seen in (p!7|), the coupling 
of the chemical potential explicitly breaks Lorentz invariance. This will play an 
important role in our analysis of the quantum equivalence. 

^One can also investigate different linear combinations of the Noether charges [M. 
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To perform the perturbative analysis first one has to determine the classical 
ground state of the system, in other words the minima of the Hamiltonian 
Ti = 'H{9, (j), iIj^po^p^^p^)^ finding the critical points and checking the Hessian's 
positive definiteness. Using the actual (perturbative) expression for r, eq. (p^), 
it is straightforward to show that in the perturbative region around r = — 1 
(which in the original model corresponds to (7 = —1 as opposed to the dual 
model where it is exactly -1) there is a two-parameter family of local minima, 
all of them being physically equivalent, given by: {pg = Pip = 0,p^ = j,6 = 
^,(p = const. ,ip = const.). Since (j) and ^p does not appear explicitly in the 
Hamiltonian, we can chose for convenience the stable classical ground state, 
common to both models, to be given hy 9 = ^, (p = tl^ = 0. We shall expand 
our fields around this background. 

The Lagrangian, C, as emerging from eq. ([l5| ) is not suitable for a direct 
perturbative computation on account of the overall 1/A factor. However this 
factor can be removed by an appropriate rescaling of the fields: 

= j, <A=5, i^ = ^, (18) 

where A = y/X as it was given in (jl^). We note here that originally the fields 
had no renormalization in either the original or in the dual model, consistent 
with the fact that the Euler angles are compact variables. The rescaling nev- 
ertheless introduces nontrivial renormalization, but this can be deduced from 
the renormalization of A. 

Deleting the tilde from the rescaled fields, we get the following result: 

^ =lmOf + (1 + rcos^{Ae))id^(P + ^5o,f+ 

h 

+ {d^ipf + 2aujf,u cos{Ae){df,(l) + -6of,)du^j} . 

According to the general prescription of Wick rotation, we define the Eu- 
clidian continuation of our model by Ce = —^nido — > ido,di — > di). Note 
that a similar continuation of h ^ ih would lead to inconsistencies (like wrong 
signs in the propagators: j_^2 instead of ) ■ 

Before starting the perturbative expansion of the trigonometric functions 
we make another field transformation (redefinition): 9 ^ ^ + 9, guaranteeing 
that the minimum we expand around is: {9 = 0,(p = 0, "0 = 0). Thus the 
Euclidian Lagrangian we use has the form: 



^ =l{i^^,9f + (1 + rsm\A9)){d^(t> + ^<5o^)'+ 

+ (0;,^)^ + 2asm{A9)iu;^,ud^,<j)duip + ^J^o!.)^^^} , 



(20) 



where r and a are given by ( p^q ) resp. (|14|), and the new parameters w and $7, 
which bear the model dependence, are: ujfj,u = —S^u, = iSou in the original 
and cj^jy = —ie^y, = —eqi, respectively in the dual model. 

Next we expand the Lagrangian C around the classical ground state: {9 = 
0,(/) = CV' = 0), with A = y/\ being the relevant coupling constant. We note 
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at this point that one could follow a different route and expand the parameters 
r, g and a in terms of the RG invariant. The motivation for this would be that 
at the end of the computation this has to be done anyway. This possibility and 
the complications that arise will be investigated elsewhere ||20t| . 

After some algebra the result is as follows (for the sake of simplicity we 
denote £ by £): 

C = C-2A-^ + C.iA-^ + £o + +CiA + C2A^ + o{A^) , (21) 

where 

C-2 = -^h^, C-i = -ih do(j) , 

A =lid^i^f + lid^^^f + lid^Of - rhiH'' + a/iO, Qd^i, , ^^^^ 
L\ = - irh O'^doip + auj^u 9d^(t)du'4' , 

£2 =]-re\d^(t>f + \rh^9^ - lahn.e'^d.^p . 

2 b b 

Notice that £_2 is a constant (i.e., it is independent of the fields), while 
C-i{h,ip) = —ihdo(j) is a total derivative, thus it can be discarded in this 
non-top ological sector of the theory. From Cq we see that (p is a massless scalar 
field, while 9 and ip are mixed, apart of the mixing the former is a massive 
scalar field with mass \/—rh, the latter is massless. The interaction of the dif- 
ferent fields is highly non-trivial, as can be seen above, and contains infinitely 
many vertices. Nevertheless, these vertices are naturally separated in the weak 
coupling regime. We emphasize that only the first two terms, C-2 and £-1, are 
common to both models, as the model dependent parameters Qj, /?, uj and O 
appear in Cj for all j > 0. 



4 The free energy 

Our goal is to define and compute the free energy density in perturbation theory. 
After setting the stage we do the explicit computations. 



4.1 Definition of the free energy 

At this point our aim is to define the free energy density perturbatively. De- 
noting the fields collectively hy ip = {9,(l),'ip), the free energy (density) reads: 

where V is the volume of the system, S[h,^p] = J(Px C{h, ip{x)), and jVip 
denotes the functional integration over the field configurations ip = {9,(j),%l)): 
jDip = J'D(j)'Dil;'D9. The role of the denominator in ( ^3|) is to insure the cor- 
rect normalization: J^{h = 0) = 0. From dimensional arguments one expects 
the following functional dependence: J-{h) = h?'^{h)^ where ^ = '^{h) is a 
dimensionless function. 
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Let us note the similarity between the free energy defined above and quan- 
tum effective action (the generator of the IPI graphs). The role of the external 
field is played by /i, that couples to a conserved charge (composite operator) 
rather than an elementary field. This similarity will play a structurally simpli- 
fying role when we discuss the renormalization of the model. 

For the perturbative expansion, in view of eqs. (21,2^, it proves useful to 
introduce 



oo 



S[h,^]= Si[h,^]A\ (24) 

i=-2 

with Si[h^(p\ = JdPx Ci{h,if{x)), where i > —2. Using this in eq. ( p^ ) we 
obtain a similar expression for the free energy: 

oo 

J'ih) = J2 :Fi{h) A\ (25) 

j=-2 

Our task will be to determine the first few terms in this expansion, in both 
models, and compare them. More precisely we determine the first six terms of 
J-{h) and check whether they are equal. 

As £_2 is independent of the fields ip, it results that S-2[h,ip] = —^h'^V, 
and exp(—S-2[h,(p]) factorizes (we will come back to this) in the functional 
integral (|2^). This way one readily obtains that the first term of ( |25| ) is 

jr^,(h) = -^h\ (26) 

Of course this is valid both in the original and dual model, implying that at 
leading order the perturbative free energy densities coincide. In addition, since 
£_i is a total derivative, it implies that S-i[h, ip] = 0, thus 

T-i{h) = 0; (27) 

again a model independent statement. 

Thus all what remains to be dealt with is the reduced action 

oo 

S[h,^]=Y,S,[h,^]A\ (28) 
and the reduced free energy 

oo 

^ih) = J2^iih)A\ J'ih) = r^2{h) + ^{h) . (29) 

Eq. (^) might also suggest that all J^i with i odd vanishes, this would comply 
with the fact that the coupling A = is just an artifact of our perturbation 
theory, as it was A that appeared in the original Lagrangian. The vanishing of 
all odd power contributions in A would indeed imply that the true coupling is 
in fact A. 
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Introducing the Z(h) = J'Dipe ^^^'"^^ auxiliary function then we can rewrite 
eq. (23) as e^'^^^'^^ = Z{h)/Z{h = 0). Moreover, if M is an operator we define 



the following 'expectation value': 

Expanding Z{h) as a power series in A we have: 
Z{h)=[l-{S,)A-{S2-]-S,^)A^- 

^ r _ (31) 



{S3 - SiS2 - ^S^)A^ + ^(A^) ]jv^e 



For simplicity we have omitted to write the functional dependence of the 
Si[h,(^]-s. Using the identity 

1 + yl + y2 ^2 + 2/3 ^3 + o(^4) = ^+"^ ^'+"3 , (32) 

where ^ ^ 

xi = yi, X2 = 2/2 - -Vi , X3 = y3-yiy2 + -yf , (33) 

we can read off the various components of the reduced free energy density 

T,ih) = ^{Si), (35) 

:F2{h) = ^[{S2-ls,') + l{S,)']. (36) 

We emphasize that in the above formulas we kept only those terms which depend 
on h. In other words we discarded the contribution of Z{h = 0), which in fact 
is a divergent quantity. This is consistent with the fact the that only the 
derivatives of the free energy are observable thermodynamical quantities. 

4.2 Propagators 

To compute the various vacuum expectation values (or correlators since we are 
in Euclidean space) determining J-i{h) we use dimensional regularization (with 
n = 2 + e). Since 

SqQi) = J (fx£o{h) = ^J (fxip\x)M{x)^{x) (37) 

where 

99* (x) = (0(x),0(x),V(x)), 

-(9^ - r/i2 ahVL^dy 
M{x) =1 -d^ I , 9" = d^^df, , 
—ahO^udu —d"^ 



(38) 
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in order to determine the propagators of the various fields we have to invert the 
matrix operator M(x). This is easily done in momentum space resulting: 

Geix) = e-^P-p^gip) , G^x) = e-^P- (p^ - rh^)gip) , 

Ge^x) = iahj^ e-^P-n.pg{p), G^^{x) = e'^f^ ^ , 

Ge^{x) = 0, G.^^{x) = 0, 

(39) 

where 

We have used commonly the notation G^p^'{x — y) = {if{x)ip'{y)), with G^ = 
G^^. Just note in passing a few simple properties: Ggg{x) is even, while Gg^{x) 
and G^^{x) are odd, and G^q{x) = —Gq^{x). It is also obvious that (j) behaves 
like a massless scalar. 

4.3 Computation of J-'o{h) 

The computation of J-Q{h) involves in fact the evaluation of a functional deter- 
minant, similarly to the the case of the quantum effective action. Evaluating 
the Gaussian integral from eq. (34), using the identity detX = e^"^ "^^^^ ^ results: 



1 r (Pp _ h^{rp^ + { aVL ■ pf) . 

2 J (27r)" ^ 



Mh) = ^ / T^e^ In (1 - " ■ (41) 



A proper way to compute this expression is to take its derivative with respect 
to h and solve the following initial value problem: 

d^„W_ ^„(0) = 0. (42) 



dh J (27r)" p4 _ /j2(^p2 _^ (^an ■ pf 

Rescaling p — > hp, the h dependence factorizes and we get: 

T(h\ = -— f "^"^ rpp^ + (gpO • p)^ 

n J (27r)" p4 _ h2^r^p2 + ^aon ■ p)^] " ^ ^ 

Above we have made it explicit that the integral is computed in terms of the 
bare quantities rg and oq, rather then the renormalized ones, r and a. Prom 
now on we are going to make this distinction clear in all subsequent formulas. 



This same expression was obtained in |17], though the initial Lagrangian 
differed from the one used here by a certain rescaling of the fields. A rescaling 
usually cannot cause major discrepancies, at least at low orders of perturbation 
theory, as this example also refiects. 



The model dependence is manifest in (42). The integral is divergent in 
two dimensions more precisely it has a first order pole in e. The analysis of 
f[% has computed the pole term and the constant term in the e expansion of 
(^) in closed form, in terms of generalized hyper-geometric functions. After 
renormalization, the two expressions - as functions of the original (respectively 
dual) coupling and parameter - were not equal, but, the difference could be 
accounted for by the scheme dependence of the two loop beta functions. Indeed 
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it was pointed out in Q that the equivalence of the original and dual (3 functions 
corresponds to a perturbative redefinition of the coupling constants in the dual 
model: 

A = A + ^(l + <7), g = g+^il + gf. (44) 

Implementing this redefinition in the expressions of the renormalized one loop 
free energy densities revealed their equality. 

Here we take a different route from the one described above. Our strategy is 
to express the parameters that bear the model dependence in terms of the RG 
invariant in both models, then set these two RG invariants equal and compare 
the results. Evaluating —r^p^ — {aQ^L-pY we get — goPi in the original model 
and Pi — gopQ respectively in the dual model, where Pfj_ = e^uPu- Due to |26| we 
have p^ = P^i and one can perform a change of variables from p to p. This way, 
changing also pi <-> PO) one can obtain formally identical expressions in the two 
cases. At this point we only remark that the role of po and pi in dimensional 
regularization is different. 

4.4 Computation of J-'i{h) 



According to ( |3q ) and (|22|) the computation of J-i{h) involves the following 

correlation functions: {9{x)di,ip{x)) and {6'^{x)do(j){x)) . Using Wick's theorem 
we obtain: 

J^i{h) = 0. (45) 

Eq. (^) has twofold meaning. It shows once again the model independence of 
the free energy, though we have already encountered explicit model dependence. 
On the other hand supports our earlier statement about the vanishing of the 
non-analytic corrections in A. 

4.5 Computation of !F2{h) 

In the case of J-2{h) we will not be able to obtain the result in a closed form, 
nevertheless what we can actually compute will suffice to achieve our goals. 
According to eq. (36) we have to compute {S2{h)), (S'i(/i)^) and (S'i(/i))^. We 



can immediately quote eq. (45) and {Si{h)) = 



4.5.1 Computation of {S2{h)) 
From eq. (^2|) we see that 



{S2{h)) = I d^x [^ro{e^{x){d,,cpf{x)) + -roh^{e\x))- 
1 



b 



(46) 



At first sight one might want to discard the terms that are not coupled to 
h. Nevertheless these terms acquire h dependence through the /i-dependent 
propagators. 

Based on Wick's theorem for the first term we have {9'^{x){d^(j))'^{x)) = 
(0^(x))(((9^0)^(x)). Since the 4> propagator is in fact a Green's function, or in 
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other words fundamental solutions of the corresponding wave equation, we can 
conclude that {{d^(j)f'{x)) = (5*^"')(0) (modulo equal time commutator terms), 
where 5^"^ (0) is the Dirac delta distribution 'evaluated at 0'. But in dimensional 
regularization the latter is set to zero, yielding no contribution. 

The next term in eq. (H) is also readily evaluated: (6'^(x)) = 3(6'2(x))2 = 
3G5)(0)^. Evaluating the third term gives (recall that the system in finite vol- 
ume: / dTx = V): 



{S2ih)) = ]fV 



n 



GeiO? . (47) 



4.5.2 Computation of {Si{hf) 

The computation of (S'i(/i)^) involves in fact the evaluation of double integral 
jdP-xdP-y {C{x)C{y)). From eq. (p2|) it turns out that (5i(/i)^) equals jiPx 
times 

rlh^diG^{x)[Ge{Qf + 2Ge{xf] - 4irohaoojf,udod^G4x)Ge{x)d^Ge4x)+ 
+ aQUj^„uJxpd^,dxG^{x)[-d^Gg^{x)dpGg^{x) - Go{x)d^dpG^{x)+ 
+ d^Ge^{Q)dpGe4,m- 

(48) 

A priori it is not clear at all how to obtain the overall factor required by 
dimensional analysis. Moreover it is also puzzling how the explicit factors of i 
disappear during the computation. As we shall see shortly it is the form of the 
'tensor' parameters Vl^, lo^i, and propagators that is responsible for the correct 
answers. 

The first term in (|48| ) results 

-2r^h^V-GeiOf6oo. (49) 
n 

The only non-trivial fact that one has to use an IR regularization of the field 
(j), the regulator's mass we denote by m. Then the first term in this parenthesis 
will be proportional to 

d^^xd^GAx)^ AW(fc) = 0. (50) 

J fc^ + 

The second term in (Hsl) results 



Aroa^h VtiVp^np J j^j^ {h+W^ k, g{ki)k2uk2p9{k2) , (51) 

while the last terms of (^) give 

.2,. "^"^2 {ki+k2)x{ki + k2)p 



ao^^^^XpV / ^7^^^)2 9{k,)9{k2)k2p ^^^^ 

[a^h^ki^ki ■ n){k2 ■ n) + k^k2u{ki - ro/i^)] . 
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The general structure of the integrals that appear (with one exception) is 
of the following form: 



The only integral that cannot be brought to this form has k2g{k2) instead of 
g{k2). If g{k) were a covariant expression in fc, then the value of the integral 
would be given completely by the index structure. More precisely covariance 
would require the result to be the sum of triple products of the Kronecker 
delta functions. The number of independent possibilities is 6!/(2!)^/3! = 15. 
Exploiting the obvious symmetries under the exchange of fJ-i ^ f^2, ^J'3 ^ M4, 
^ fJ-G, and the less obvious exchange (/^3,^4) (^Si/^e); which can be seen 
by a change of integration variables, the tensor structure reduces significantly 
to only four independent terms. Denoting Sf^^^■ by (ij), we get for (^) 




(12) (34) (56) /i + (12)[(35)(46) + (36)(45)]/2 + (34)[(15)(26) + (16)(25)]/3 + 
+ (56)[(13)(24) + (14)(23)] h + {(13)[(25)(46) + (26)(45)] + (14)[(25)(36)+ 
+ (26)(35)] + (15)[(24)(36) + (23)(46)] + (16)[(23)(45) + (24)(35)]} I4. 



Using the standard methodology one can compute the unknowns Ii through I4, 
in a straightforward manner. Unfortunately g{k) is not covariant with respect 
to the full SO{n). As we pointed out, covariance was already broken at the 
level of the Lagrangian. If we analytically continue po to uq dimensions and 
pi to ni dimensions, with uq -\- ni = n, then instead of the full SO{n) group 
we get S'O(no) x 5*0(711). In other words, among the po-s and pi-s standard 
covariance arguments remain valid, and the computation sketched above makes 
sense. Thus if all the indices are of po or pi type, then we can reliably compute 
the integrals. 

In the expression of the two-loop free energy (^, [5^ ) the free indices we 
have dealt above are contracted with different tensor structures. In the original 
model and uj^u involve only delta functions, and this way what we said 
applies. On the other hand for the dual model and iOfj^u are epsilon tensors, 
mixing the indices, and the arguments presented break down. The same is true 
for the sum of the terms involved, that we eventually want to compute. 

There is an independent argument that shows that even if we were able to 
do the integrals, the result would not be reliable, due to the behavior of the 
e tensor in dimensional regularization. More precisely, even if the integrals in 
(H) and dH) were covariant, based on the definition of u;^,^ and on (U), we 
could conclude that in the dual model (5i(/i)^), and as a result J^2{h), contains 
terms proportional to the product of two e tensors, with uncontracted indices. 
The broken covariance by the external field /i, invalidates the above argument, 
but is highly likely that it complicates matters, rather than simplifies, and as a 
result we would still end up with e^j^ea/j terms. As it is well known, there is no 
consistent way defining such on object in dimensional regularization |]2^ . 

At this point the computation of (5*1 (/i)^) is hopeless, driving the same state- 
ment about J-2{h). Nevertheless we can do something less ambitious. Following 




(53) 




(54) 
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g, in order to compare the results, at the very end of the computation we 
want to express them in terms of the two RG invariants, which are finahy set 



equal. In other words, based on (|T^) and (16) we trade the renormalized quan- 
tities g and r for the RG invariant M, implicitly meaning that by this we have 
also set the RG invariants equal. Since it is the bare that appears in (^) 
and (|5^, and also in (42), let's investigate more closely what happens to this 
term during renormalization and expansion in terms of the RG invariant. 

At the level of the bare quantities we aim to have oq = \/l + go- Based on 
Q this is certainly true at two loop level in the original model, while in the 
dual one this is more subtle. As it was shown in Q naively the relation oq = 

+ go cannot be maintained at two loop level in the dual model. Nevertheless 
taking into account the redefinition of the dual model's coupling and parameter, 
eg . (]4^) , the above relation can be maintained at one loop order. 

As ao = \/l + go is doubtlessly valid in the leading order, the renormal- 



ization (|5|) amounts to ao = yl + Zg(X,g)g. But (^) shows that Zg{X,g) = 
1 + o(A), thus to leading order we have = -^/l + g, with g the renormalized 



coupling. On the other hand at this point we can use (|llD and conclude that, 
after renormalization and expanding in the RG invariant, ao becomes propor- 
tional to Vx. 

The good news is that all the terms we were unable to compute (^ij) and 
(^) are proportional to a^, hence are of order A^. Since we have not even 
attempted to compute the J^A^h) term that is of the same order, we neglect 
them for the time being. Having in mind the insertion of the RG expressions 
for the parameters, we simply get: 

^2ih) = -h^'^GeiOf . (55) 

Naturally this result is to be interpreted as modulo terms that will be of higher 
order after renormalization and expansion in the RG invariant. 

5 Renormalization 

The results from the previous section are divergent, and we used dimensional 
regularization to compute them. As advertised we make use of the similarities 
between the free energy and the quantum effective action to discuss the renor- 
malization of the former, modeled by the renormalization of the latter (see e.g. 
j21[|). In this section we follow a more or less naive renormalization procedure, 
and compute the renormalization group improved perturbative two-loop free 
energy. It will be the role of the next section to tight the loose end, and prove 
that what we did is indeed correct. 

Since we are at second order of perturbation theory, we have both first and 
second order poles in dimensional regularization, as can be seen in eq. (^). 
Since the free energy is a physical quantity, it has to be well defined after 
renormalization. The recipe of this section is to use the renormalization of 
the deformed principal a model (^),(|^, and the renormalization of its dual, to 
renormalize the free energy. The above renormalizations were performed using 



the geometric method of [23, 
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This procedure can be immediately objected since the free energy is com- 
puted in a theory that has additional terms in the Lagrangian (^), compared 
to the deformed principal a model (Q) and its dual (^. A priori there is no 
reason to expect that the wavefunction and coupling constant renormalization 
functions are the same, though in fact they are. In the following we present 
a simple argument in favor of the above statement. We shall give a complete 
proof in section 5.5 below. 

Coupling the external field amounts to the appearance of terms proportional 
to h. Setting /i = we get the original models (deformed principal a model 
and its dual). Thus the wavefunction and coupling constant renormalization 
functions can at most differ from the ones of the original models by terms pro- 
portional to h. But /i is a dimensionful quantity, having mass dimension +1 (it 
is a super-renormalizable coupling or relevant perturbation). On the other hand 
our theories have no other dimensionful quantities, and the wavefunction and 
coupling constant renormalization functions are dimensionless. We conclude 
that h cannot appear in the latter ones, proving the assertion. 

In addition since h couples to a conserved charge, it is not renormalized. 
The only issue that remains will be to deal with the renormalization of the 
operators that couple to /i, viewed as composite operators. We postpone this 
to section 5.5. 

In order to cancel the second order poles of the regularized free energy (^) 
we have to go beyond the computation of Let us review what we know at 
this stage about the renormalization of the deformed principal a model and its 
dual. In (|l|) we have introduced a generalized Lagrangian, and argued in the 
last paragraph of the section that it encompasses the renormalization properties 
of both models. Thus we can translate the renormalization properties of the 
original models to those of the generalized Lagrangian. The Euler angles {(j)fi,tp) 
are compact and have no wavefunction renormalization. The renormalization 
of r, A and a in the unified model follows from those of g (resp. a) and A in the 
original models. 

In ( [l8|) we rescaled the fields {(l),9,ip), this way in (|2^ ) they have the cor- 
responding nontrivial wavefunction renormalization. Thus Z\(X,M), the cou- 
pling constant renormalization of A, is the central object for their renormaliza- 
tion. (Recall that M is the renormalization group invariant that appeared in 
(|lO|)). Zx{X, M) was computed at two loop order in perturbation theory in Q. 
The interest in was restricted to the single pole terms. On the other hand 
we are constrained to deal with the second order poles too. First we compute 
the residue of the second order pole in Z\{X,M). 

5.1 A at two loop 

The goal is to determine the terms yx{X,M) and yx{X,M) in the following 
expansion: 

Za(A, M) = 1 + - yxiX, M) + 1 yx{X, M) + o(l). (56) 

The first term was basically determined in Q. All we need is to use @) for 
Zx{X,g) and ( |T^ for g, in terms of the renormalization group invariant M, to 
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obtain: 



(57) 



Our task of computing y\{\,M) is highly simphfied by the special properties 
of the non-linear a models. Following |22] it was shown in [23| that the gen- 



eralized renormalization theory |24] of non-linear a models lead to generalized 
renormalization group equations, that allow one to determine the residues of 
the higher order poles in a given coupling constant renormalization function 
like Z\{\,M), without extra diagrammatic computations. 

More precisely it is shown that in a theory (like ours) with a single di- 
mensionless coupling constant A, with mass scale parameter defined in n 
dimensions (dimensionally regularized and minimally subtracted), having an 
expansion of the bare coupling Aq in terms of the renormalized coupling A of 
the form 



An 



A + 



oo 



-(n-2)-;' 

the pole residues ay{X) satisfy the recursive pole equations: 



d 



(l_A|-)ai(A)|^a,(A). 



(58) 



(59) 



We can apply this directly to compute yx{X,M). From the expressions 
following (^) we see that in our case ai(A) = A yx{X, M), with y\{X, M) given in 
(p7|). Using ( [5g| ) for = 1 we can determine the leading term in the expansion 
of a2(A). A simple calculation shows that 02 (A) = t^A^ -|- . . ., implying 
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yA(A,M) = ^A2 + ... . 



(60) 



5.2 One loop free energy 

By now it is a computation on the back of an envelope to obtain the renor- 



malized one loop free energy density. We have to use (26) for the leading term 
(with A as the bare coupling Aq) and ( ^3|) for the one loop regularized result. 
The integral itself, as it was pointed out in |17|, is hard to deal with exactly. 
Nevertheless, with the parameters traded for the RG invariants, and keeping 
only the leading contribution in A (as explained in the last paragraph of Section 
4), and renormalizing the expression using (^), ( |56| ) and (57), we obtain: 



l~loop 



'P{h) = - — 



1 

2A 



1 

8^ 



In ( - 



7 - 1 - ln(47r) 



(61) 



Some of the terms that are higher order in A, and are not yet displayed, will be 
used for the computation of the J^2- 

We note two things here. Firstly, the one loop equivalence might be argued 



to be not surprising based on the experience gained in |1C] and the general 



one- loop beta function result of [11|. Secondly, the above expression correctly 
reproduces the known one-loop free energy of the 0(3) sigma model, that arises 
in the ai = limit. 
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As opposed to the corresponding computation in ||T^, the procedure that 
leads to (^) does not require the sphtting of the momentum integration in the 
0-th direction and the rest. The reason is simply that due to the expansion 
in the RG invariant, the order A term becomes totally covariant. The non- 
covariance is shifted to the next order. 



5.3 Two loop free energy 



For the computation of the two loop free energy we need all the results devel- 
oped. The o(A) corrections arising from the leading bare term (|2^ ) are obvious 
in the light of @, (|56|), ( |57| ) and (|60|) . The o(A) contributions from the bare 
next to leading term of the one loop free energy ( ^3|) can be computed as an 
expansion in the RG invariant, using (14) and ([T^). The results are: 



Xh 



2«1 



2 r 



1 

47re 



1 

8^ 



(7-ln(47r) +2 ln(/i) - 1) 



(62) 



where Ai equals Jn in the original model, and eou^Ou in the dual one. Although, 
as argued above, the product of two epsilon tensors with uncontracted indices 
is ambiguous in dimensional regularization, the contraction of one index gives a 



meaningful expression. Following |26|, we assume that e^ae^a has a consistent 
continuation, namely: 



'-UfJ, 



-'fJ.U • 



(63) 



These two expressions might differ depending on the regularization schemes 
chosen. As we have pointed out in connection with (54), we can a priori continue 
the 0th direction into no dimensions and the 1st direction into ni dimensions, 
provided uq + ni = n. We checked that the final result is consistent with 
dimensional analysis in either case. While we know |17] that the choice uq = 1 



is a consistent scheme with the continuation of the e tensor, we cannot claim 
the same about the general case. 

As we will see in a moment the only discrepancy in the two-loop free energies 
of the two models comes from (p2|) . This way we have several choices of different 



schemes to see the quantum equivalence. With the notation following (62) we 
can chose in both models the scheme with Jqo = ^ii- The second choice is two 
different schemes in the two models, but related to each other by: 6qq = 6ii. 
Besides these naive choices we have a highly non-trivial one, with the identical 
choice in both models: (5oo = 1 and 6ii = n — 1. Unlike for the previously 
mentioned ones, we know the consistency of this scheme, and we are going 
to work with it in the rest of the paper. As we remarked in connection with 
(^), it was shown in [17| that in this case the difference in (|6^ ) is accountable 
for a perturbative redefinition of the coupling constant A in the dual model. 
Naturally this redefinition has to proceed the expansion in the RG invariant. 
Unlike in |17], where this gave correction to T^^''°°^{h), due to the expansion 
in the RG invariant we get contribution to T'^^^°"^{h). The redefinition of the 
coupling constant does not effect the genuinely higher order terms as it can be 



seen from (44) 
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Thus, as we just noted, we are going to use the last scheme in what follows: 
(5oo = 1 and 6ii = n — 1, and for definiteness we consider the case of the original 
model first. In the light of what we just said, for the dual model, the results 
(U), (|65| ) below, must be amended by ^^'^^ to account for the perturbative 

redefinition of A, eq. (^4|). However once this redefinition is taken into account, 
the two loop results below are identical in the two models. The result for the 
other schemes will differ from the results to be presented by terms that are 
some number times — A/i^^. 



Using the bare result ( |55D and summing it with the corresponding contri- 
butions described, we obtain a cancelation of both the first and second order 
poles, and obtain a finite expression. The computation is somewhat tedious, 
but straightforward, the result is simply: 



^ ^ 2A Svr 



In ( - j + 7 - 1 - ln(47r) 



167r2 



In 



+ 7 



ln(47r) 



In 



+ 7 - 1 - ln(47r) 



(64) 

It simplifies a bit if we use instead of minimal subtraction (MS scheme) the MS 
scheme: In// — > In/x + (7 — ln(47r))/2. Then the free energy reads: 



2— loop 



167r2 



/l2 










Svr 










{ 


In - 






V/V 



1 



In ( - 



(65) 



5.4 Improvement of the perturbation theory 

We can take advantage of the asymptotic freedom of our theory and calculate 
the RG improvement of the perturbative result. This gives the asymptotic 
expansion of the free energy for large values of the external fields. 

Physical quantities depend on the renormalized coupling A, the renormalized 
parameter 5, and the dimensionful scale parameter (or subtraction point) //, in 
such a way that the action of the renormalization group (RG) operator 



2^ = /.| + /3a(A,,)|^ + /?,(A,,)|, 



(66) 



vanishes on them. As the free energy density takes the form J-{h) = 
—h'^'if{X,g,fi,h) and the external field h is not renormalized, the function ^ 
is renormalization invariant = 0. Since we are interested in the behaviour 
of the free energy density for large values of h we write h = Hqc^ where ho is 
fixed and let x — > 00. Standard RG considerations then give that 



J^{h) = -h''^{X{x),g{x),f,,ho), 



(67) 



where A(x) and g{x) are the running coupling and parameter. In our final 
result for the two loop free energy density (^) the parameter g is eliminated 
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in favor of the two loop RG invariant. As a consequence, up to this order, 
^' = '^{X{x),fj,, ho) depends only on A(x), which is a solution of 



^\{x) = /3a(A(x), ai), \{x = 0) = A, 
where the beta function has the following expansion: 



(68) 



/3a(A) 



-6oA^ - 6iA^ - ftaA^ + 



(69) 



with (see eq. (12)) 
" 27r' 



6, = ifi 



47r V vr 



1 

4^ 



2 



27r 



(70) 



Thus we can go on with the RG analysis as if we had only one coupling constant 
A(x). The expression of the (RG invariant) A parameter is 



A 

or more conveniently 



e >'o^X ''0 e' 



1 + 



bi hi 



A + o(A2 



(71) 



In 



A 



1 / 
M ~ hi 



4 In A + c + 



h± 



A + o(A2 



(72) 



/i UQy 

where c is a constant of integration. We define A = A^ by choosing e'^ = 

^) ''o as it simplifies some of the forthcoming expressions. 

As the next step of the RG analysis we note, that using expression (|72[), it is 
customary to define an effective coupling a{h) by the following transcendental 
equation |28|: 



hn In 



h 1 6i 
; = — h Ina. 



^MS " ' ^0 



(73) 



The important property of this effective coupling is that it depends on the 
physical quantity s = In ( - J^ ) , moreover one can express the running coupling 

A in terms of the effective coupling a perturbatively: 



A = a (1 + iia + ^2"^ + 



(74) 



where 



6 = 6o ( In — - - c 



^ c ^2,^1. hi h2 



(75) 



To use the o(A) part of the free energy density effectively we need the large s 
asymptotic expansion of the effective coupling up to o(s~^): 



27r / ^Ins B Cln^s Dins 
a = — IH + — + ^ + 

S \ S S S'^ s 



E 



2 + s2 



(76) 



where the coefficients A, . . . E can be obtained from ([7^). 



22 



Now improving ( |65| ) by the RG, i.e. using eq. (74-7^) in ( |65| , [67| ), results, 
when the dust settles, the following formula for the asymptotic (large s) be- 
haviour of the two loop free energy: 



(/.)=--(. + (l-|)ln.-- + (l-|)^- 



1 

+ - 

s 



1 - - + 



S2 



(77) 



Unfortunately 62 has not yet been computed in the literature, nevertheless we 
can say a lot about it. Due to the expansion in the RG invariant, it will have 
terms coming from the lower order beta function coefficients, like in (^), and 
the genuine three loop coefficient evaluated at g = —1. For simplicity we denote 
this last term h^^^ = b^^~^°°P\g = —1). But this is exactly the three loop beta 
function coefficients of the 0(3) cr-model: b^^^ = b^^^\ This beta function has 
been computed in [29|, and in our convention for the coupling constant b^^^ 



takes the form b^^^ = 3^^. Carrying one step further the computation of (^) 
results: 

^-8^("N^?)+"''-T^(f-3p) + 3^. (78) 
Accordingly, as a final result we have: 

47rV 2 2 2 s s 4 s^/ 

5.5 Renormalization of the composite operators 

In this section we give a solid foundation to the results of the previous renor- 
malization procedure. More precisely we consider the renormalization of the 
composite operators that arise via the coupling of the chemical potential. Our 
attitude is similar to the standard procedure of mass renormalization in QCD: 
initially the mass of the light quarks is set to zero, then the bilinear mass opera- 
tor il^ip is added, and its effect is accounted by its renormalization as a composite 
operator. We followed the same ideology so far, neglecting the fact that the 
renormalization of the coupled theory is different from that of the uncoupled. 
It is now that we remedy this. 



Let us reconsider the computation of the free energy (23). In (24) we started 
with the bare Lagrangian. Instead we have to use the renormalized one, with 
counterterms coming from the wavefunction renormalization of the fields, cou- 
pling constant renormalization, and the renormalization of the composite oper- 
ators, all these coming with a natural grading: 

00 

S[h, if] + AS[h, J2 (^^[^' + ^^dh, V^]) A\ (80) 

i=-2 

with the renormalized quantities in the right hand side. ASi[h,{p] commonly 
denotes all the counterterms of order i. The original Lagrangian that has to be 
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employed reads: 

^° = lid^^f + \idt.(t>? + lid^O? + Aau;^, 05^09,^ + A^^r9\df,^f . (81) 

AS'_2 amounts simply to the multiplicative renormalization we considered in 
(^). Since multiplied by a factor of 1/A, it gives rise to a term that is inde- 
pendent of A (already used to renormalize the one loop free energy), one that 
is o(A) (also used), and higher order terms: 

AS-2 = ^S% + AS^_2X + o{X^). (82) 

AS"-! = as already S-i = 0. A priori ASq contains terms from the from the 
wavefunction renormalization of the fields, but these are independent of h, and 
are canceled by the denominator in (p3[). Thus ASq = J^o + 6qi O^, where 
OO ^ _ri/i26)2 g^^f^ Qi ^ ahn^ediyt/;. The renormalization of 0° (and the 
similar operators that appear in the Lagrangian proportionally to h) has two 
contributions: one from the renormalization of h (that is zero as discussed in the 
previous section) and one from the renormalization of the composite operator 
9^. We expect that 6qo = 1 + o(A^), and this expectation can be confirmed by 
a short explicit computation. 

In order to compute 5qo we need a Green's function involving O^. For 
simplicity let's consider the one-point function of 6'^: {9'^{x)). As a renormal- 
ization prescription we normalize the one-point function according to the tree 
level value, and determine 5q2 from the condition of preserving the above nor- 
malization. At o(A) there are two diagrams that contribute coming from the 
o{A) and o{A^) vertices. The latter one has a value proportional to 

j cTy Gix - yf d^G{^) = J d^y G{x - yf 5{G) = (83) 

in dimensional regularization. The same result is obtained if the diagram is 
evaluated in momentum space, where the masslessness of the fields requires 
additional IR regularization (as exploited already in a previous section). 

The other diagram is readily proportional to \a?. The remaining integral 
can be computed but we don't need the result for what follows, because once 
again, as we express this contribution in terms of the RG invariants we shall 
have a dependence proportional to A^. Thus we have 5o = ^ + °(^^) least. 
As a consistency check we quote that the same result is obtained by considering 
for example the {9'^{x)9{y)9{y')) Green's function. 

During the computation we employed the ^h^O"^ operator as a mass term 
for 9, though it is a composite operator as discussed above. The motivation for 
this can be given as follows: assume that is a perturbation, and expand 

it perturbatively with the rest of the terms in (|3l|). The difference is that 
^h?9'^ is independent of A, and as such any term in its power series expansion 
is of the same order, and has to be summed. In other words, any term in 
the perturbative expansion of J- will be multiplied by the full expansion of the 
exponential of ^h'^9'^, that can be resummed. The resummation on the other 
hand is equivalent with the corresponding mass term h for the 9 field. From 
the point of view of the original h = theory this is a non-perturbative result. 
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Next we consider the renormalization of the second operator O and the 
one appearing at the next level o{A): = O'^d^cf). These operators will have 
their composite operator renormalization functions: 5qi and 5q2. Since the 
undeformed {h = 0) theory has only interactions proportional to A or to a^/X, 
and the latter ones must appear at least twice for a finite contribution, eq. (^), 
we conclude that the renormalization effects due to these vertices are at least 
of order A and resp. Xo?: 6qi = I + a^X + . . . resp. = 1 + z'^X + . . .. In 
the spirit outlined above we have to introduce the new terms z^O^ and z^O^ 
into the action ( ^0|) as terms contributing to AS[h,(f], and account for their 
contribution. As these new terms are proportional to A, they contribute to 
AS2[h,(f]. Thus we expect new contribution to J^2- Based on (|3T|) these are 
proportional to (z^O^) resp. (z^O^). But we already know from ( |45| ) that these 
are zero. It is easy to see that all the other combinations of these two operators 
with the rest of the operators will give higher order contributions to the free 
energy. 

It is obvious that the product of z^O^ (resp. z'^O^) with (resp. O^) 
will give non-zero contribution to J^^, but as far as J^2 is concerned there is no 
deviation from the results we obtained. 

At the next level (i.e. at o(^^)) we have the composite operators = 6^ 
and = O^di^Tp. Consider first = 9^. In order to compute 6q3 we fix 
the normalization of the correlation function {9'^{x)) to the tree level value: 
{6*{x)) = 3(5(0)^. We compute 5q3 by computing the first corrections to the 
correlation function, and insisting as above that the above normalization re- 
mains valid. Using the available interaction terms in (^) we see that using the 
vertices Odfj^cpd^il^ resp. 0^(9^(/))^ give the first nontrivial contribution propor- 
tional to A^ resp. A^, while their combination is of order A^. Thus we conclude 
that 6q3 = l-|-o(A^). In the case of after simple considerations it is similarly 
obtained that 6q4, = 1 + o(A^), and we are done. 

From the above analysis we conclude that the more or less naive renormal- 
ization procedure employed in the first place is completely adequate, and the 
additional composite operator renormalizations do not disturb the results as far 
as the J^2 is concerned. 

6 Conclusions 

In this paper we have computed the change in the ground state energy density 
of a deformed principal SU{2) sigma model, and one of its T-duals, in the 
presence of an external field. The computations has been carried out at two 
loop order in perturbation theory. Perfect agreement has been found in the 
following sense: there were renormalization schemes in the two models that 
yielded the same expressions for the two loop free energy densities. 

We defined the free energy using the conserved charge corresponding to a 
special symmetry of the model. Care had to be made at the choice of the 
charge, since we wanted it to exist as a genuine symmetry in both the original 
and dual model, while it is known that the charge corresponding to the isometry 
used for the duality becomes topological. We performed the computation for 
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the simplest choice, but other possibiUties are also investigated at least in the 



one loop order [27|. We computed the relevant diagrams and performed 
the renormalization using dimensional regularization. Due to asymptotic free- 
dom, we could improve the perturbative result using the renormalization group 
equations. 

This result strengthens the confidence in the findings of Q , that were mainly 
based on beta function computations. More precisely this way we have given 
much stronger evidence that in the case of constant goo, at least in the case 
of the deformed principal SU{2) sigma model, the naive Buscher formula that 
relates the original model and its T-dual, gives a true quantum equivalence. As 
it was pointed out and exemplified in [^], the constancy of goo is no guarantee 
for two loop quantum equivalence. Based on this one might expect that a 
finer test would detect discrepancy even in the case where the beta function 
arguments show no sign of it. Though the free energy test that we performed 
does not prove the lack of discrepancies, it hints to the non-existence of these 
at least at two loop order. It would be interesting to analyze the pertinent 
SL{3) example of 110] from the free energy point of view, though, at this point 



it seems to be too big of a computational challenge. 
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